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ABSTRACT 

It is currently believed that the "standard" accretion disk theory under-predicts the observed X-ray luminosity 
from Soft X-ray Transients (SXT) in quiescence by as much as 4 to 6 orders of magnitude. This failure of the 
standard model is considered to be an important argument for the existence of the alternative mode of accretion 
- Advection Dominated Accretion Flows (ADAF) in astrophysics, since these flows allow a much higher level 
of X-ray emission in quiescence, in agreement with the observations. Here we point out that, in stark contrast 
to steady-state standard disks, such disks in quiescence (being non-steady) produce most of the X-ray emission 
very far from the last stable orbit. Taking this into account, these disks can accommodate the observed X-ray 
luminosities of SXTs rather naturally. Our theory predicts that Fe Ka lines from standard accretion disks in 
quiescence should be narrow even though the cold disk goes all the way down to the last stable orbit. 

Subject headings: accretion, accretion disks — binariesxlose — novae, cataclysmic variables — line: formation 



1. INTRODUCTION 

Soft X-ray transients are mass transfer binaries which peri- 
odically undergo outbursts in which their luminosity increase 
by several orders of magnitude (e.g., see reviews by Tanaka & 
Lewin 1995; van Paradijs & McClintok 1995; Wheeler 1996; 
Cannizzo 1998). The compact object may be either a neu- 
tron star (NS) or a black hole (BH). SXTs are particularly in- 
teresting given the ongoing debate about the structure of the 
inner accretion disk. In quiescence, the luminosity of SXTs 
seems to be much higher (e.g., Verbunt 1995; Lasota 1996, Yi 
et al. 1996; Lasota, Narayan & Yi 1996a; Narayan, McClin- 
tock & Yi 1996) than allowed by the standard model (Shakura 
& Sunyaev 1973; expected to be modified for SXTs due to 
hydrogen-ionization instability - e.g., Mineshige & Wheeler 
1989 [MW89]). Lasota, Narayan & co-workers argued that 
because the "standard" theory predicts that the inner accre- 
tion disk is essentially empty of mass in quiescence and ac- 
cretes at a very small accretion rate, M m , it cannot repro- 
duce the observed X-ray luminosity £ x ,obs: typical numbers 
are 0.1M in c 2 ~ (1(T 4 - 1CT 6 ) x L Xj0bs . They conclude that 
SXTs in quiescence clearly invalidate the standard model. 

However, in this Letter, we argue that this argument is appli- 
cable only to "bare" standard accretion disks - those without an 
overlying hot corona. We find that quiescent standard disks with 
a corona should produce most of their X-ray emission at very 
large radii precisely due to the fact that the inner disk is almost 
empty of mass. At large radii, the standard theory of hydrogen- 
ionization instability allows much larger accretion rates, and 
the observed X-ray luminosities can be easily reproduced if the 
corona reprocesses ~ 10% of the disk energy budget. We also 
make a testable prediction - Fe Ka lines from quiescent ac- 
cretion disks should be narrow (even though these disks extend 
down to the last stable orbit or the surface of the compact star). 

2. HYDROGEN IONIZATION INSTABILITY: BASIC FACTS 



Because the argument of Lasota (1996) against the standard 
model in SXTs is qualitative, we will address it with a qual- 
itative model as well. The hydrogen ionization instability of 
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accretion disks has been studied for more than two decades for 
a range of central objects - White Dwarfs (WD), NS, and galac- 
tic and extra galactic BHs (e.g., Hoshi 1979; Smak 1982; 1984; 
MW89; Mineshige & Shields 1990 [MS90]; Cannizzo 1993; 
Siemiginowska, Czerny & Kostyunin 1996 [SCK96]). Can- 
nizzo (1998, §2 & 3) gives a short overview of the development 
of the theory and an easily accessible and enlightening discus- 
sion of the background physics. Here we will only sketch the 
basic principles of the instability to the extent needed to explain 
our main point. The behavior of the unstable disk depends on 
the so-called "S-curve" that relates the local disk effective tem- 
perature and the column density (see Figures in Smak 1982, 
MS90, SCK96; and Fig. 31 in Frank, King & Raine 1992). 
The curves are different for different radii, but luckily only 
slightly. The steady state Shakura-Sunyaev solution predicts 
that the disk effective temperature, T c s(R), changes as i?~ 3 / 4 
with radius. In general, this predicts that for small R, hydro- 
gen is completely ionized, then there should be a range of radii 
where hydrogen is partially ionized, and finally, at large R, hy- 
drogen should recombine completely. 

However, it turns out that the partially ionized solution is 
thermally unstable. In terms of the S-curve, this region is be- 
tween points C and E (see Fig. 31 in Frank et al. 1992). Nu- 
merical simulations show that the unstable disks oscillate be- 
tween two quasi-stable states. The accretion rate is not con- 
stant with time and radius, but the time-averaged accretion rate 
equals that at which the matter is supplied at the outer edge of 
the disk, Mo. One of these states - the region between points 
E and D on the S-curve - is the outburst state in which hy- 
drogen is completely ionized. The other thermally stable state 
- the one between points B and C on the S-curve in Fig. 31 
of Frank et al. 1992 - is the quiescent state when hydrogen 
is neutral. To remain in this state, the mid-plane gas tempera- 
ture should be kept below ~ 10 4 Kelvin, which in turn requires 
T eS (R) < T C (R) < few x 10 3 (depending on the black hole 
mass, M, and viscosity prescription). 

Time-dependent (full disk) calculations of many authors 
show that when the unstable region of the disk returns to quies- 
cence, it settles around point B on the S-curve, which is a little 
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below the maximum quiescent temperature and surface density 
(point C). The accretion rate corresponding to point B, M B (R), 
is lower than M , so that the matter starts accumulating, most 
quickly at the outermost part of the unstable region. At a future 
time, this piling up of the mass will lead to M(R) > Mq at 
some radius (equivalently, T c g(R) > Tq{R)), and the outburst 
will be triggered at that radius (see Cannizzo 1998, §3). 

Fortunately, for the present analysis, these details are mostly 
irrelevant. What is important for us here is that both local disk 
calculations and full disk time-dependent ones show that the 
gas effective temperature in quiescence is between T B (R) and 
T C (R); that T B (R) is typically smaller than T C (R) by a fac- 
tor of 2 or less, and that both of these functions are extremely 
weakly dependent on R. In particular, the scaling law given 
by Smak (1984), MS90 and SCK96 is T C (R) oc R~ q , where 
q is between 0.08 and 0.12. Full disk calculations give tem- 
perature profiles slowly increasing and then decreasing with R 
(e.g., Fig. 13 and 16 in MS90; Fig. 10 in SCK96), which 
is probably due to the fact that the outer disk fills up quicker 
and hence lies closer to point C than the inner disk does. 
There are also observations of the disk effective temperature 
in quiescence for eclipsing dwarf novae, systems to which the 
hydrogen-ionization instability was applied first. These obser- 
vations show that T c s(R) is either a constant or a slowly de- 
creasing function of R (Fig. 1 1 in Wood et al. 1986; Fig. 14 in 
Wood et al. 1989; and Figs. 5 & 6 in Bobinger et al. 1997). 

While the exact run of T c g(R) with radius depends on de- 
tails of the model and the object under consideration, all of the 
above observational and theoretical papers show that T e g(R) ~ 
const in quiescence within a factor no larger than 2. This is 
quite remarkable given that the radii of the unstable region 
can span some two or more orders of magnitude: the steady- 
state Shakura-Sunyaev solution with T e fi(R) oc i?~ 3 / 4 would 
predict a change in T c g by about a factor of 30. Hence, be- 
low we will assume that T e g(R) ~ Tu = const ;$ Tc ~ 
few x 10 3 K in quiescence. Physically, the near isothermal- 
ity of the disk in quiescence is due to the fact that the ioniza- 
tion state of hydrogen is a very sensitive function of tempera- 
ture, and large changes in radius cause only a slight difference 
in T. For example, the opacity is an extremely strong func- 
tion of midplane disk temperature T^: in the temperature range 
6000 < T d < 8000 K, the opacity oc T h d where b is as large as 
10-20 (see Cannizzo & Wheeler 1984). 

Of course, the effective temperature is approximately con- 
stant only inside the unstable disk region. Let us call the cor- 
responding radius Ru- For larger radii, the disk is always non- 
ionized and so it follows the steady state Shakura-Sunyaev so- 
lution, i.e., T c s(R) oc R~ 3 ' 4 . Summarizing, 



T eS (R) 



Tu, 

T SS (R) 



if R < Ru 
otherwise , 



(1) 



where Tss(R) is the Shakura-Sunyaev effective temperature of 
the disk corresponding to the given radius, accretion rate M 
and the black hole mass. The thermal disk flux is simply given 
by F d (R)=a B T c %(R). 

What are the implications of this for the accretion rate 
through the disk? At R > R s = 2GM/c 2 , F d (R) = 
{3/8tv)GMM(R)/R 3 , which can be solved for M(R): 



M(R) 



M (R/Ru) 3 , 
M 



for R< R v 
for R> Ru 



That is, the accretion rate scales as M(R) oc R 3 within the un- 
stable radius. Physically, this comes about from the fact that it 
is most difficult to keep the hydrogen recombined for smaller 
radii because the gravitational and viscous energy release scale 
as R~ 3 , and hence the disk has to be most strongly depleted for 
the smallest R to remain quiescent (see also Cannizzo 1998). 

3. X-RAYS DO NOT ALWAYS COME FROM THE INNER DISK 

A frequent implicit or explicit assumption in the common 
logic about the X-ray emitting region(s) is the following. Since 
X-rays are the "hottest" part of the overall disk spectrum, and 
since the Shakura & Sunyaev (1973) model predicts that the 
disk is hottest in its innermost region, then that has to be the 
region where the X-rays are produced. Seemingly innocent, 
this suggestion is actually appropriate for the ADAF-type mod- 
els (i.e., the inner-ADAF plus outer cold disk models), but is 
actually entirely foreign for standard disks with coronae in qui- 
escence. 

To show this, let us first assume without any justification that 
the corona above the disk exists at all radii, and that the fraction 
of power transferred into the corona, /, is independent of radius 
(e.g., Svensson & Zdziarski 1994). That is, F X (R) = fF d (R). 
Next, we define the X-ray luminosity per decade of radius, as 
L X {R) = R 2 F X (R), which obviously scales according to 



L X (R) 



fLu(R/Ruf 
fL u (R/R u y 



for R < R v 
for R> Ru 



(3) 



where Lu — F d (Ru)R v . Therefore, most of the viscous en- 
ergy liberation and X-ray production for a quiescent disk occurs 
at R ~ Ru 3> 3Rs, if / = const. The inner disk in the stan- 
dard model is simply very dim and for many purposes may be 
unobservable at any wavelength in quiescence. 

Is it appropriate to assume / ~ const? To date, detailed 
MHD simulations of accretion disks have not yet addressed this 
question, because global MHD disk simulations are forbiddenly 
numerically extensive, and hence calculations are limited to a 
small region in radius (see, e.g., Miller and Stone 2000). In the 
absence of these, we will make a simple analytical argument. If 
one defines the a-parameter through v = ac s H, where c s and 
H are the sound speed and the pressure scale height, respec- 
tively, then the total energy generated by the disk per unit sur- 
face area is 2F d (R) = (9/8) {GM/R 3 )vY, = (9/8) ac s P toU 
where E is the vertical mass column density, and P tot is the to- 
tal disk mid-plane pressure. The flux of magnetic energy, F mag , 
delivered by buoyancy to the disk surface and then presumably 
transformed into X-rays by reconnection is i^mag ~ vi,(P mSLg ), 
where Vb ^ c s is the average buoyant velocity, and {P ma g) is 
the average magnetic field pressure (~ energy density). The 
viscosity parameter a was defined by Shakura & Sunyaev 

(p ) 

(1973) to be a ~ a t + x " ag/ , where the first term, a t is 

■Hot 

the turbulent viscosity coefficient. It is now believed that the 
magnetic viscosity is significantly more important than the tur- 
bulent one (e.g., Hawley, Gammie, & Balbus 1995; 1996), and 
that a ~ (Pmag)/-Ptot, and hence 



mag 



v b aP to t ^ ac s P t( 



F d 



(4) 



(2) 



In other words, this equation shows that / is indeed expected to 
be a constant with respect to R. 

Further, if corona is powered by magnetic flares (e.g., Galeev, 
Rosner & Vaiana 1979; Haardt, Maraschi & Ghisselini 1994; 
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Beloborodov 1999; Miller & Stone 2000), then it is clear that 
the resulting X-ray spectra have little to do with the tempera- 
ture of the underlying disk. The spectra depend primarily on 
the Thomson depth of the hot plasma and its temperature. The 
temperature within a flare is controlled by the ratio of the X-ray 
flux to the disk flux, F x /F d , (e.g., Stern et al. 1995; Pouta- 
nen & Svensson 1996), plus the exact (unknown) geometry of 
the reconnection site. Therefore, it is entirely possible that a 
magnetic flare that occurs at R ~ 1000 Rs will have an X- 
ray spectrum similar to, or even harder than that of a flare at 
R ~ 10i?5, especially for the quiescent disks. 

To summarize this discussion, there is nothing special about 
the innermost region of the disk for quiescent accretion disks. 
If there is a corona above the disk, then most of the emission 
at any wavelength will be produced at R ~ Ru [Note that in 
a steady Shakura-Sunyaev disk, most of the viscous energy re- 
lease happens in the inner disk, and this is why it is safe to 
assume that X-rays come from very small radii for such disks]. 
Since there are several parameters that will determine the re- 
sulting X-ray spectra (e.g., reconnection rate, height and ge- 
ometry of the X-ray source, Thomson depth of the X-ray emit- 
ting plasma, etc.), it is difficult to calculate the spectra model- 
independently. Nevertheless, we note that variations of these 
parameters are sure to yield a broad range of spectral indices 
(see, e.g., Stern et al. 1995; Poutanen & Svensson 1996), and 
hence it may be entirely possible to reproduce the observed 
spectra of SXTs. 

4. X-RAY "PROBLEM" FOR STANDARD DISKS IN QUIESCENT SXTS 

Currently, it is believed that there exists a clear-cut "evi- 
dence" for a major failure of the standard disk instability model 
for the quiescent state accretion disks around dwarf novae, neu- 
tron stars and galactic black holes (SXTs). The popularity of 
the ADAF model itself in part rests on this fact because it is 
claimed to be the only model that can explain self-consistently 
the level of X-rays in quiescence for the systems mentioned 
above. Let us for a moment accept the suggestion that X-rays 
do emanate from the innermost region and follow the arguments 
of Lasota (1996), as well as many of other authors (e.g., Lasota 
et al. 1996a; Yi et al. 1996). The requirement for the disk ma- 
terial to be below the instability point C on the S-curve limits 
the accretion rate to be (Lasota 1996): 

M m (R) ~ 2.76 x 10V r?- u Mf °- 37 «~ a79 g s" 1 , (5) 

where Mi is the mass of the black hole in units of 10 Mq, t% 
is the recurrence time in units of 10 9 s, and = i?/10 7 cm. 
Numerical models (e.g., MW89) satisfy this requirement. This 
condition (eq. 5) is similar to the simpler one given by equation 
(g) with a scaling M m (R) oc R 3n instead of R 3 . Using a 10 % 
efficiency for conversion of mass accretion rate into X-rays, La- 
sota (1996) noted that two SXTs, A0620-00 and V404 Cyg, 
have been detected at the accretion rates ~ 10 10 and ~ 3 x 10 12 
g s _1 (see Verbunt 1995). This estimate is 4-6 orders of mag- 
nitude larger than that given by equation (^). 

However, as we argued in §|l most of X-ray emission comes 
from R ~ Ru ^S> 3Rs for standard disks in quiescence. Ap- 
proximately, we have 

tijj ZKu 

where M(Rjj) ~ M . Since in quiescence 



M(Ru) ~ M m (Ru), we estimate 

L x < 3.7 x 10 23 ftg 1 rf "M?- 63 a~ - 79 erg s" 1 . (7) 

Clearly, by allowing rj = i?[//10 7 cm to be much greater than 
unity, we can obtain much higher L x than the estimate based 
on 10% efficiency and r ~few. 

Let us now use theoretical and observational constraints on 
the parameters in equation (g) to determine L x more accurately 
for A0620-00. From observations, we know that L x ~ O.lLbb, 
where Lbb — 10 32 erg s _1 is the blackbody optical luminosity 
of the disk as inferred by McClintock et al. (1995). This means 
that we need to choose / ~ 0.1. Further, analyses of Ha and 
Hp spectroscopic data (Marsh et al. 1994) indicate that the ra- 
dius where the Balmer line luminosity peaks is ~ 5 x 10 10 . 
This should roughly equal Ru since that is where most of the 
luminosity is produced in our model. Therefore, 

° 3 °(aTs) (^fV»f. (S) 

where ol\ = a/0.1. This value of L x is in a very good 
agreement with the observed L x ~ 10 31 ergs _1 . A further 
independent check of the self-consistency of this picture can 
be made by computing the temperature Tjj, from <JbT v — 

(1 - f)(3/&n)GMM /R$ T . For M ~ 10 16 g s" 1 (see Mc- 
Clintock et al. 1995), we find 

/ M \ 1/4 

r -= 3 - 2xi ° 3 Ud ■ (9) 

This temperature estimate agrees quite well with the expected 
temperature at the lower kink of the S-curve (see Figs, in Smak 
1982, MS90, SCK96). 

5. DISCUSSION 

In this paper we re-examined the argument given by Lasota 
(1996), Lasota et al. (1996a) and Yi et al. (1996) with which 
they rule out the standard accretion disks extending down to 
the last stable orbit for quiescent states of SXTs. The argu- 
ment is based on the assumption that all X-rays are emitted 
very close to the last stable orbit. Because the accretion rate 
through this region is so low in quiescence, the standard theory 
appears to under-predict the observed X-ray luminosities by or- 
ders of magnitude. However, we have shown that this argument 
is only applicable to bare standard accretion disks, i.e., those 
that have no hot X-ray producing corona. Such disks are ruled 
out already due to the fact that they do not produce hard X-rays 
at all, and hence there is no reason to invoke them as a model 
for SXTs in the first place (and not only for SXTs but for all X- 
ray emitting accreting sources). We have shown that if a corona 
exists at all radii, then most of X-ray emission will come from 
large radii in quiescence. At these radii, the maximum allowed 
accretion rate - the maximum one at which hydrogen is neutral 
- is much larger. Accordingly, much larger X-ray luminosi- 
ties may be sustained by such "non-bare" disks, in qualitative 
agreement with the X-ray observations of quiescent SXTs. 

Further progress in testing this and the ADAF models for 
SXTs can only be made by a detailed comparison of theoretical 
and observed spectra. While there is a large body of literature 
on the spectra of the ADAF model for SXTs (e.g., see the re- 
cent review by Lasota 2001) and none for the model proposed 
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here, both of these face serious theoretical challenges and hence 
neither model can be preferred (in our opinion). In particular, 
standard disks with magnetic flares is perhaps the most viable 
model for luminous accreting black holes, both AGN and GB- 
HCs (e.g., Galeev et al. 1979; Beloborodov 1999; Nayakshin & 
Dove 2001; Done & Nayakshin 2001). However, it is not clear 
at all how spectra of the high luminosity magnetic flares (e.g., 
Poutanen & Svensson 1996) should be modified when applied 
to SXTs. Ironically, the uncertainties in the theory may be great 
enough to preclude both positive or negative statements about 
the "standard" model for SXTs, if these require detailed spec- 
tra. 

The ADAF model, on the other hand, has fundamental diffi- 
culties in the modeling of the cold-disk ADAF transition, which 
is therefore often done in an ad-hoc way. Another major un- 
certainty of the theory is the accretion rate through the disk 
as a function of radius (M (R)). The earlier papers assumed 
M(R) = const. However, this would only be the case if the 
ADAF solution joined the cold disk at R t £ Ru, where M(R) 
is indeed constant (see eq. ^). Yet R t is treated as a free param- 
eter of the model with M being independent of it (e.g., Lasota 
et al. 1996b; Yi et al. 1996). From our analysis above, we be- 
lieve that if one chooses R t < Rjj, the accretion rate through 
the ADAF region should be smaller by the factor (R t /Ru) 3 
(see also Wheeler 1996, on a similar point). Since the ADAF 
luminosity scales as M 2 for small accretion rates, this factor 
translates into the sixth power of R t /Ru in L x . In addition, 
the issue is further complicated by the more recent discovery of 
the importance of convection and winds in the ADAF solution 
(e.g., Blandford & Begelman 1999), neglected in the earlier pa- 
pers. Taken together, this implies that at the present, the X-ray 
spectra of the ADAF model are no more well determined than 
the spectra of magnetic flares. 

One obvious consequence of our model is that because most 
of the X-ray power is radiated away at R ~ Ru, Fe Ka lines 
of quiescent SXTs will have to be narrow even though the disk 
extends down to the last stable orbit. While narrow-ish lines are 
also predicted by the ADAF model, we believe that a detailed 



modeling and observations of these lines are still the key to dis- 
tinguish the two models. The ADAF model produces the lines 
mostly due to collisional processes, and the line is dominated by 
H- and He-like components (e.g., Narayan & Raymond 1999) 
with the line centroid energy of ~ 6.9 and 6.7 keV, respectively. 
In the case of a cold disk plus corona model, depending on the 
parameters of X-ray spectra we can expect the line centroid en- 
ergy to be either at ~ 6.4 or 6.7 keV or be a combination of 
these two (e.g., Nayakshin & Kallman 2001). Highly ionized 
collisionally excited Fe Ka lines may also be expected based on 
the analogy with Fe Ka lines from solar flares (e.g., see Pike et 
al. 1996). The ADAF model is however unlikely to produce a 
neutral-like Fe Ka line, because most of X-ray emission in the 
ADAF is expected at small radii which yields little illuminating 
flux at large radii where the cold disk is located. 

Finally, we would like to stress the following general "the- 
orem". The absence of emission from the inner disk (whether 
this emission is the Fe Ka line emission, the X-ray continuum 
or the thermal disk emission), is not an evidence for the ab- 
sence of the standard disk there. The particular case of SXTs 
discussed here proves this point. Because hydrogen is neu- 
tral in quiescence, the inner disk is expected to be very dim 
at all wavelengths and so it may be difficult to discern in the 
data. This point is quite general as it may also be relevant 
for other than hydrogen ionization disk instabilities. For ex- 
ample, Nayakshin, Rappaport & Melia (2000) have devised a 
specific model for the radiation-pressure driven instabilities in 
the galactic micro-quasar GRS 1915+105. In this model, the in- 
ner disk is also very much dimmer during the "quiescence" than 
expected in the steady-state Shakura-Sunyaev solution (see Fig. 
13 in Nayakshin et al. 2000). Therefore, one must not use the 
steady-state Shakura-Sunyaev theory for unsteady disks - e.g., 
low luminosity systems, such as SXTs and Low Luminosity 
AGN (since one may suspect that the reason these objects are 
dim is that they are in a quiescent state), and those verifiably 
unsteady such as GRS 1 9 1 5+ 1 05 . 

The authors acknowledge many useful discussions with John 
Cannizzo, Chris Done and Demos Kazanas. 
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